A signature for the Luttinger liquid phase in alternating Heisenberg spin-1/2 chains 
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We have studied the zero and low-temperature behavior of anisotropic alternating 
antiferromagnetic-ferromagnetic Heisenberg spin-1/2 chains in a transverse magnetic field. Using 
the analytical spinless fermion approach, the thermodynamic behavior of the model has been stud- 
ied. We have introduced a new order parameter to distinguish the gapless Luttinger liquid phase 
from the other gapped phases. The exact diagonalization Lanczos results are also used to compare 
with the spinless fermion ones. We have found a double peak structure in the specific heat curves 
for the region between the two quantum critical points. Using the numerical full diagonalization 
results, the existence of the double peak structure is confirmed. 

PACS numbers: 75.10.Jm, 75.10.Pq 



I. INTRODUCTION 

The study of continuous phase transitions has been one 
of the most fertile branches of theoretical physics in the 
last decades. Each phase can usually be characterized by 
an order parameter. Often, the choice of an order pa- 
rameter is obvious. However, in some cases finding an 
appropriate order parameter is complicated. In partic- 
ular, many magnetic systems experience the Luttinger 
liquid (LL) phase for certain values of some non-thermal 
control parameter. Luttinger liquid is the paradigm for 
the description of interacting one-dimensional (ID) quan- 
tum systems^* 2 -. The correlation functions decay as power 
laws and the ground state has a quasi-long range order. 
In spite of all manifestations of this phase, the magnetic 
ordering in the LL phase is still unclear, i.e, all suggested 
local magnetic order parameters turn out to be zero, hav- 
ing no effect in determining the structure of the LL phase. 
In this work, we will concentrate on studying the proba- 
bly nonlocal magnetic ordering of the LL phase. 

One dimensional quantum spin systems such as anti- 
ferromagnetic spin-1/2 chains 3 -, spin ladder systems^, or 
bond alternating spin-1/2 AF-F chains are good candi- 
dates for studying the LL phase. When these systems are 
placed in an external magnetic field, they can be mapped 
onto a ID system of interacting spinless fcrmionsi££i£&2. 
The filling of a fermionic band can thus be continuously 
tuned, making these systems suitable for probing the LL 
physic o 1 ' 3 ' 10 ^ 1 . In the case of isotropic bond alternating 
AF-F spin-1/2 chains the whole band can be covered. In 
the absence of an external magnetic field (h = 0), the 
model is mapped onto a nonlinear sigma model with a 
Airs topological angled. This model is always gapful and 
can be regarded as a Haldanc gapped spin-1 chaini 3 -. The 
gap decreases by increasing the magnetic field and goes 
to zero at the lower critical field h Cl —. At h Cl , the system 
enters the gapless phase and the fermionic band starts to 
be filled. This process continues until the upper criti- 
cal field h C2 is reached. Increasing the field beyond h C2 
reopens the gap and the band is then completely filled. 

Recently, It was demonstrated that CuBr^CsH^N^ 



is a unique system for controlling and probing the physics 
of LLi£. This sample is a spin ladder system with 
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14.6 Ti&. Bond alternat- 



ing AF-F spin-1/2 chains that allow experimental ac- 
cess to the whole fermionic band are not known. How- 
ever the (CH3)2NH2CuCl3 system has been considered 
to be a suitable realization of bond alternating AF-F 
spin chain. Linkcd-cluster calculations and bulk mea- 
surements show that DMACUCI3 is also a realization 
of the spin-1/2 alternating AF-F chain with nearly the 
same strength of antifcrromagnetic and ferromagnetic 
couplings^!. Other experimental samples of the AF-F 
alternating spi n-1/2 chain compounds have also been re- 
ported in Refs.fHIHiimiH. Althou g h [t is not known 
whether these systems can experience the LL phase, the- 
oretically, an isotropic bond alternating AF-F spin chain, 
in its ground state, can enter the LL phase with quasi- 
long range order. In this paper, we look for an order 
parameter to describe the LL phase of the bond alter- 
nating AF-F spin-1/2 chain. The Hamiltonian of this 
model with anisotropic ferromagnetic coupling is given 
by 

N/2 

N/2 N 

i=i 3=1 

where S", (a = x, y, z) are spin-1/2 operators on the j-th 
site. Jf and Jaf denote the ferromagnetic and antifcr- 
romagnetic couplings, respectively. A is the anisotropy 
parameter and h is a uniform magnetic field. 

The ground-state propertie a 23 ' 24 i 25 ' 26 i 27 ' 28 i 29 and low- 
lying excitations 3 ^ of this model have been thoroughly 
investigated by numerical tools and variational schemes. 
In particular, the string order parameter which was orig- 
inally defined for the spin-1 Heisenberg chains 3 ^ has been 
generalized to this system. The ground state has long- 
range string order, which is characteristic of the Haldanc 
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phase. Hida has shown that the Haldane phase of the 
AF-F alternating chain is stable against any strength of 
randomness^. The ground state phase diagram of the 
AF-F alternating chain in a longitudinal (a = y) mag- 
netic field is studied using numerical diagonalization and 
finite-size scaling based on conformal field theory^. It is 
shown that the magnetic state is gapless and described 
by the LL phase. The model represented by the Hamilto- 
nian in Eq. [1] which includes a transverse magnetic field 
has been studied recently using the numerical Lanczos 
methodil The main attention in this study was focused 
on the investigation of field-induced effects in the ground 
state phase diagram, present when the antiferromagnetic 
coupling Jaf dominates {Jaf > Jf)- The system has 
two critical fields and the energy gap in the intermedi- 
ate region depends on the anisotropic parameter A. For 
A / 1, the intermediate state is gapful and the ground 
state of the model has stripe-antiferromagnetic order—. 

In this paper we consider again an anisotropic AF-F 
chain in a transverse magnetic field. Using the numerical 
exact diagonalization method and the analytical spinlcss 
fermion approach, we investigate the zero temperature 
and thermodynamic behavior of the model. We introduce 
a new mean field order parameter which can distinguish 
the LL phase from the other gapped phases (Fig. 
In the specific heat curves versus temperature, a double 
peak structure appears in the intermediate region of the 
magnetic fields h Cl < h < h C2 . 

The outline of the paper is as follows. In section II 
we discuss the zero-temperature ground state phase dia- 
gram of the model. In section III we present the results 
of the spinless fermion approach and the numerical full 
diagonalization results on the low-temperature behavior 
of the model. Finally we conclude and summarize our 
results in section IV. 



II. ZERO-TEMPERATURE BEHAVIOR 

A. RESULTS FROM THE NUMERICAL 
LANCZOS METHOD 

In this section wc briefly discuss the model (1) in the 
limiting case of the strong AF coupling Jaf 3> Jf- In 
this limit the model can be mapped onto an effective 
spin-chain Hamiltonian^. At Jaf 3> Jf, the system 
behaves as a nearly independent block of pairs^. Indeed 
an individual block may be in a singlet or a triplet state 
with the corresponding energies given by 



Ei,-i 



Jaf 



Jai 



3Jai 



For h < Jaf , one component of the triplet becomes closer 
to the singlet ground state such that for a strong enough 
magnetic field we have a situation when the singlet and 
S z = 1 component of the triplet create a new effective 
spin r = 1/2 system. On the new singlet-triplet sub- 
space, the original Hamiltonian becomes the Hamiltonian 



of a fully anisotropic XYZ spin- 1/2 chain in an effective 
magnetic fielcU^ 

J'=l 



N/2 

i=i 



(2) 



where h eff = h - Jaf + Jf/4. At A = 1, the effective 
problem reduces to the theory of the XXZ chain with a 
fixed antiferromagnetic anisotropy of 1/2 in a magnetic 
field^. The gapped Haldane phase at h e f f < h*f f = — 4f 
for the AF-F alternating chain corresponds to the neg- 
atively saturated magnetization phase for the effective 
spin chain, whereas the massless LL phase of the AF-F 
alternating chain corresponds to the finite magnetization 
phase of the effective spin-1/2 chain. The critical field 
\\ e JJ = Ie. where the AF-F alternating chain is totally 
magnetized, corresponds to the fully magnetized phase 
of the effective spin chain. 

Away from the isotropic point A = 1 the effective 
Hamiltonian @ describes the fully anisotropic ferromag- 
netic XYZ chain in a magnetic field that is directed per- 
pendicular to the easy axis. In this case, it is founcU^ 
that a gapped stripe-antiferromagnetic phase exists for 
the intermediate values of the transverse magnetic field 
h ci < h < h C2 . 

To recognize the different phases induced by the trans- 
verse magnetic field in the ground-state phase diagram, 
we implemented the modified Lanczos algorithm of finite- 
size chains (N = 12, 16, 20, 24)H with J AF = 1, J F = 1/2 
and different values of the anisotropy parameter A. The 
energies of a few lowest eigenstates were obtained for the 
chains with periodic boundary conditions. In Fig. [1] wc 
have plotted the results of these calculations for different 
values of the anisotropy parameter A = 1.0, 0.5, and 
chain length N = 20. As it is clearly seen from this fig- 
ure, in the case of zero magnetic field the spectrum of the 
model is gapped. For h ^= the gap decreases linearly 
with h and vanishes at the critical field h Cl (A). In the 
isotropic case A = 1.0, the spectrum remains gapless for 
h ci = 0.78 ± 0.01 < h < h C2 = 1.0 and becomes once 
again gapped for h > h C2 . But, in the anisotropic case 
A = 0.5, the excitation spectrum is gapfull except at two 
critical field values h Cl = 0.82±0.01 and h C2 = 0.94±0.01. 
In the intermediate region h Cl < h < h C2 the spin gap 
which appears at h > h Cl , first increases vs external field 
and after passing a maximum decreases to vanish at h C2 . 

In conclusion, at T — 0, two quantum phase transi- 
tions in the ground-state phase diagram of the model 
have been identified with increasing transverse magnetic 
fielcU^. The first transition corresponds to the transition 
from the gapped Haldane phase to the gapless LL phase 
(or gapped stripe-antiferromagnetic for the anisotropic 
case) . The other one is the transition from the gapless LL 
phase (or gapped stripe-antiferromagnetic) to the fully 
polarized phase. 
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spinlcss fcrmions; 




FIG. 1: The excitation gap of a spin-1/2 AF-F chain versus 
the uniform magnetic field h, for different anisotropy param- 
eters A = 1.0, 0.5 and chain size N = 20. 



FERMIONIZATION 



The behavior of the gap and the quantum phase space 
of the model lead us to investigate the thermodynamic 
properties of the model. In this respect, we imple- 
ment the Jordan- Wigner transformation to fcrmionizc 
the model. Because of two types of coupling constant, 
we introduce two kinds of spinlcss fermion through the 
following Jordan- Wigner transformations: 



9+ 

°2n-l 



ale 



2 n 



bit 



°2n-l 



a'a r , 



$2n — blpn 



Using the above transformations, the thermodynamic be- 
havior of the isotropic (A = 1) Hcisenberg AF-F spin- 
1/2 chains has been studied in the absence of a mag- 
netic fields. By the above transformations, the AF-F 
spin chain is mapped onto a ID system of interacting 
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N/2 



N/2 

£( 



n=l 
1 - A 



(°AF f. 



Jf Jaf 
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+ h.c) + ^ ^ J AF a \iCLnb\ l b'n 



N/2 

E J Fblb 
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n=l 



(3) 



Treating the Hamiltonian Hf in the mean field approx- 
imation, the interacting fermionic system reduces to 
a ID system of new non- interacting dynamical quasi- 
particles. Many mean field order parameters (auxiliary 
fields) might be considered. Many of them are irrelevant, 
i.e., they do not have a stable mean field solution, while 
some are relevant. In our system we have introduced 
the magnetization, ferromagnetic and antifcrromagnctic 
dimers as mean field order parameters; 



(a) n a n ) = d a , (blfin) = db, 
(al +1 b n ) = P F , (b\a n ) = Paf- 



(4) 



Utilizing the above order parameters, the mean field 
Hamiltonian is given by; 

Hrf = e o+^2 iz a i a k + -y&fc&k + lka\bk + h.c (5) 



where, 

A = 
B = 

Ik = 
E = 



(Jaf - J F ){d b - 1/2) + h, 
(J AF - J F )(d a - 1/2) + h, 

JAF(l/2-P AF )e lk/2 



N 

-[Jaf(\Paf\ 



Jf(P f - ~ 
d a d b + 1/4) - h 



-Jf(\Pf\ 2 - d a db + 1/4) - ^(1 - A 2 ) 2 



-ik/2 



(6) 



In the above, (. . .) represents thermal averaging over 
the Hartree-Fock eigen-states. d a and db are related 
to the magnetization, and Pp and Paf are ferromag- 
netic and antiferromagnetic exchange order parameters, 
respectively. It is clear that Paf(Pf) m the Hilbcrt space 
is the AF(F)-dimer order parameter, i.e., 



( S 2n S 2n+l) ~ ( a n+l b n) - Pf, 

(^2 n -i S 2n) = ( b n a n) = Paf- 
Using the following unitary transformations 



(7) 



a k 
bk 



u k a k + v k e l k f3 k , 
u k p k , 



v k e ltik a k 
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the diagonalizcd Hamiltonian is given by 



OLk 



k <j=± 



B-A 



V = 



i = 



(8) 
(9) 



B + A 



Equation© clearly shows that the effect of the 
anisotropy appears in the dispersion relations. The dis- 
persion relations of low-lying excitation read as, 

4 = ±[(Ja F PaF - JfPfY + JAFJFPAFpFk 2 }? + h, 

where P AF =1/2- ReP AF and P F = ±±^ - ReP F . 
Using the above order parameters, the thermodynamic 
functions such as the internal energy and the specific heat 
are expressed as follows: 



u = ^ + EE 

G 



k ct=± 



dU 



(10) 



where, nf, = \e 



is the fermion distribution func- 
tion (choose = !)• 

In order to obtain the thermodynamic behavior of the 
system, we need to know the whole temperature behav- 
ior of the order parameters for different values of h. The 
order parameters satisfy a set of self consistent equa- 
tions and it should be solved the equations ([4]), self- 
consistcntly. 

In Fig[2]-(a), we have plotted the ferromagnetic dimer 
order parameter of both anisotropic (A = 0.5) and 
isotropic (A = 1) AF-F Heisenberg spin-1/2 chains ver- 
sus h. The curves have been plotted close to the zero 
temperature. As we mentioned, at T = for different 
values of h, our ID system experiences three phases. In 
the Haldanc phase (h < h Cl ), the quantum fluctuations 
are strong enough to suppress the ferromagnetic order of 
the system and the AF(F)-dimer parameters Paf(Pf) 
should be close to the classical value —0.5(0)^. In this 
region the energy spectrum is gapped and the gap is de- 
creased by increasing h. The first quantum phase transi- 
tion occurs at h Cl . At the same time the fcrmionic band 
starts to be filled where the excited state sticks to the 
ground state and the gap of the system is closed. For the 
first region there is no specific difference between the two 
anisotropic and isotropic cases, i.e., anisotropy shows a 
marked behavior for the case of h Cl < h < h C2 . 

For the case of A = 1 and h = 0, the system is isotropic 
and has SU(2) rotational symmetry. In the intermediate 
region this system is gapless and correlations decay in a 
power law2£. In this region, the Luttinger liquid phase, 
the magnetic field suppresses quantum fluctuations and 
induces quasi long-range order in the ID system. Al- 
though the power law behavior of the correlations and 
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FIG. 2: a) The spinless fermion results of the real part of the 
ferromagnetic exchange order parameter at temperature T — 
0.05 versus h, Jaf = 1 and J F = 0.5, b) Exact diagonalization 
Lanczos results of the real part of the P F for a chain size 
TV = 20. 



other properties of the LL phase are manifest, the lack of 
an order parameter to show the behavior of the system in 
this phase is still felt. As it is obviously seen from FigJ^l 
the F-dimer order parameter has a considerable value in 
this region and is sensitive to h. The value of this order 
parameter falls sharply at the second critical field where 
the fermionic band is completely filled and the gap of the 
system is reopened. 

We have also plotted in Fig[2]-(b) the F-dimer order pa- 
rameter of the isotropic AF-F Heisenberg spin-1/2 chain 
by using the exact diagonalization Lanczos results. The 
numerical Lanczos method is implemented as a refer- 
ence approach to see the accuracy of our spinless fermion 
method. As it is clearly observed from FigJ2r(b) the re- 
sults of Lanczos method, confirms that the defined mean 
field order parameter has considerable value for the in- 
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tcrmcdiatc values of field. 

We have also plotted in Fig[3] the F-dimer order pa- 
rameter of anisotropic F-AF spin- 1/2 chain versus trans- 
verse magnetic field. Regarding the anisotropy in fer- 
romagnetic coupling, the anisotropic system is mapped 
to a ID XYZ model with U(l) symmetry. In the in- 
termediate region (h Cl < h < h C2 ) the system has the 
magnetic long range order and spins are aligned stripc- 
antifcrromagnctically. This phase is gapfull and the F- 
dimer order has no considerable value in the this phase. 
Moreover, there is a big discrepancy between the values 
of the F-dimers of the isotropic and anisotropic cases, i.e. 
in the LL gapless phase the behavior of the F-dimer order 
parameter (RePp) is different from those of the gapped 
phases. Therefore, we can conclude that the above F- 
dimer order parameter can distinguish between the gap- 
less LL phase and the other gapped phases. 

This means that experimenters might look for such be- 
havior in RePp at finite temperatures. I.e., experimental 
results on this quantity may be used to distinguish be- 
tween the isotropic chains and anisotropic ones (gapped 
stripe-antifcrromagnctic and gapless LL phases). 
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FIG. 3: The spinless fermion results of the real part of the 
ferromagnetic exchange order parameter at temperature Tf 
versus h, for Jaf = 1 and Jf = 0.5 . 



III. THERMAL BEHAVIOR OF THE MODEL 

In this section we will study the thermal behavior of 
the mean field order parameters of the Hcisenberg spin- 
1/2 AF-F chains in both anisotropic and isotropic cases. 
The order parameters also have different behavior for dif- 
ferent values of h. 

For h < h Cl , in the Haldanc phase, an increase in tem- 
perature increases the thermal fluctuations and Haldane 
ordering is suppressed. The solution of the Eqs. 2] shows 
that the AF-dimer order parameter (RcPaf) decreases 
with temperature. 

However, Jp couplings supply an interaction between 
the unit cells and try to establish F-dimer order (RePp) 
even at moderate temperatures (scaled with Jf), i-£- 
(RePp) increases with temperature up to Tp (see Fig. 
[3]). Above Tp, fluctuations are large enough to decrease 
RePp. The locations of the maximum value of (RePp) 
with respect to h are shown in Fig[3] As it is obviously 
seen from this plot, there is a peak at h very close to the 
first critical field. The location of this peak is a good 
candidate to find the first critical point. 

By increasing the magnetic field, the value of Tp (lo- 
cation of the maximum value of RePp) decreases. It is 
found that Tp goes to a minimum (Figf?]) when h in- 
creases up to h Cl . Increasing the magnetic field further 
and for h ci ~ 0.76 < h < h C2 ~ 1, Tp always maintains a 
minimum value and the value of RePp decreases mono- 
tonically with increasing T. Increasing the magnetic field 
even further, for h > h C2 , Tp also increases. 

The SF calculations also show that the behavior of Tp 
is similar to the energy gap (see FigfT]) 

The discovery of gapless or gapped excitations have led 
to the investigation of the thermodynamic properties of 
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FIG. 4: The spinless fermion results for Tf versus h, with 
Jaf = 1 and Jf = 0.5 . 



the model. One of the most important thermodynamic 
functions is the specific heat. In this respect, we study 
the temperature-dependence of the specific heat of the 
model in different quantum regimes. We have calculated 
the specific heat of both isotropic and anisotropic cases, 
using two methods, the analytical SF and the numeri- 
cal full diagonalization methods. In Fig. [5K(b) we have 
plotted the results of the SF method in the form of the 
specific heat C/N versus T with A = 1 and A = 0.5 for 
different values of h. As it is clearly seen from this fig- 
ure, in the Haldanc and paramagnetic regions (h < h ci 
and h > h C2 , respectively) , there is only one peak on the 
specific heat curve. However, in the intermediate region 
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FIG. 5: The spinless fermion results of the specific heat of an 
alternating AF-F spin- 1/2 chain versus T for different values 
of the magnetic field h with Jaf = 1 and Jf = 0.5. (a) A = 1 
(b) A = 0.5. 



FIG. 6: Full diagonalization results of the specific heat of an 
alternating AF-F spin- 1/2 chain versus T for different values 
of the magnetic field h with Jaf = 1 and Jf = 0.5. (a) A = 1 
(b) A = 0.5. 



(h Cl < h < h C2 ) a double peak is observed. At h = 
there is a Schottky like peak in the specific heat which 
is expected from the Haldane phase. By increasing h, 
the peak becomes wider and goes to lower temperatures. 
Increasing h further, causes a shoulder to appear on the 
right hand side of the curve. For h > h ci , increasing h, 
an additional peak appears which is the signal for the 
paramagnetic phase. Because of the two F and AF in- 
teractions [Jp and Jaf), there are two kinds of quasi- 
particles in the system. These quasi-particles have two 
different dynamics. These dynamics bring about two en- 
ergy scales in the system. The energy scales affect the 
behavior of the response functions. It is obviously seen 
from the plots of the specific heat C/N that the energy 
scales create a double peak in the specific heat at two 
different temperatures. As a general statement, two en- 



ergy scales can produce a double peak structure in the 
specific heat. 

We have also plotted in Figs. [B]a(b) the specific heat 
of the model by using the full diagonalization method. 
We have computed all the eigenvalues of the energies for 
different values of the transverse magnetic field h and 
anisotropy parameter A. Therefore, using these eigen- 
cnergies, we have computed the specific heat as a function 
of the temperature T. The specific heat curves have been 
plotted for Jaf = 1 and Jp = 0.5 and the anisotropy 
parameter a) A = 1, b)A = 0.5. As it is clearly seen 
from this figures, in the Haldane and paramagnetic re- 
gions (h < h Cl and h > h C2 , respectively), there is only 
one peak in the specific heat curve. To see the qualitative 
agreement between spinless fermion results and Full diag- 
onalization ones, we have plotted the specific heat versus 
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FIG. 7: Specific heat of the AF-F Heisenberg spin-1/2 chain 
vs. T for both values of A and h = 0.9 using the spinless 
fermion (SF) and full diagonalization (FD) approaches. 



T for an intermediate value of h (for example h = 0.9). 
As it is seen from Fig|7]tlie mean field results are in agree- 
ment with the others, qualitatively. It is interesting that 
the numerical results confirm the existence of the double 
peak in the intermediate region (h Cl < h < h C2 ). 



alyze the model. The first notable point is introducing 
a new mean field order parameter which can distinguish 
between a gapless LL phase and the gapped phases. This 
order parameter in the spin language of Hilbert space is 
the F-dimcr order parameter. In the isotropic case, the 
F-dimer order parameter has a considerable value in the 
LL region. We have shown that there is a big discrep- 
ancy between the values of the F-dimers of the isotropic 
and anisotropic cases, i.e., in the LL gapless phase the 
behavior of the F-dimer order parameter (i?ePp) is dif- 
ferent from that of the gapped phases. Therefore, we 
have concluded that the F-dimer order parameter can 
distinguish the gapless LL phase from the other gapped 
phases. 

The second notable point is found from the specific 
heat. We have obtained a double peak structure in the 
specific heat curves vs temperature. 

There are some questions behind the defined order pa- 
rameter which arc still controversial and we can inves- 
tigate them in future. The most remarkable questions 
referred to the topological order and spontaneous sym- 
metry breaking. I.e. in the LL phase, where the real part 
of the F-dimer is not zero, how can describe topological 
order and what kind of spontaneous symmetry breaking 
is occurred. 
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